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The  equations  of  motion  for  pressure  and  displacement  fields  in  a  waveguide  have  been  used  to 
derive  an  energy-conserving,  one-way  coupled  mode  propagation  model.  This  model  has  three 
important  properties:  First,  since  it  is  based  on  the  equations  of  motion,  rather  than  the  wave 
equation,  instead  of  two  coupling  matrices,  it  only  contains  one  coupling  matrix.  Second,  the 
resulting  coupling  matrix  is  anti-symmetric,  which  implies  that  the  energy  among  modes  is 
conserved.  Third,  the  coupling  matrix  can  be  computed  using  the  local  modes  and  their  depth 
derivatives.  The  model  has  been  applied  to  two  range-dependent  cases:  Propagation  in  a  wedge, 
where  range  dependence  is  due  to  variations  in  water  depth  and  propagation  through  internal  waves, 
where  range  dependence  is  due  to  variations  in  water  sound  speed.  In  both  eases  the  solutions  arc 
compared  with  those  obtained  from  the  parabolic  equation  (PE)  method  ©  2002  Acoustical 
Society  of  America.  [DOI:  10.1121/1.1419088] 

PACS  numbers:  43.30  Bp  [SAC-B] 


I.  INTRODUCTION 

The  coupled  mode  theory  commonly  used  in  acoustics 
was  originally  derived  by  Pierce1  and  Milder2  from  the  wave 
equation  for  the  velocity  potential.  In  this  formulation  the 
field  in  a  range-dependent  waveguide  is  expanded  in  terms 
local  modes  with  range-dependent  coefficients  (mode  ampli¬ 
tudes).  The  application  of  the  continuity  of  pressure  and  the 
vertical  component  of  particle  velocity  allows  a  partial  sepa¬ 
ration  of  the  depth  and  range  variables  and  yields  a  system  of 
second  order  coupled  differential  equations  for  the  mode  am¬ 
plitudes.  However,  as  is  pointed  out  by  Rutherford  and 
Hawker,2  while  the  boundary  condition  of  continuity  of  ver¬ 
tical  component  of  particle  velocity  is  correct  for  horizontal 
boundaries,  when  applied  to  problems  with  nonhorizontal 
boundaries,  this  boundary  condition  is  only  an  approxima¬ 
tion  to  the  correct  boundary  condition  of  the  continuity  of  the 
normal  component  of  particle  velocity.  Rutherford  and 
Hawker  showed  that  one  consequence  of  this  approximation 
is  nonconservation  of  energy  They  used  the  WKBJ  method 
to  obtain  a  solution  which  satisfies  both  the  proper  boundary 
condition  and  conserved  energy  to  first  order  in  the  slope  of 
the  nonhorizontal  boundaries  and  interfaces.  This  problem 
was  also  addressed  by  Fawcett,4  who  derived  a  system  of 
coupled  mode  equations  which  satisfies  the  correct  boundary 
conditions.  However,  in  addition  to  the  two  coupling  matri¬ 
ces,  which  is  typical  of  all  coupled  mode  theories  derived 
from  the  wave  equation,  the  equations  derived  by  Fawcett 
also  contain  two  other  so-called  interface  matrices.  These 
matrices  require  a  knowledge  of  the  range  derivatives  of  the 
local  modes,  which  can  only  be  computed  approximately. 
The  inaccuracy  resulting  from  this  along  with  the  complexity 
involved  in  solving  the  system  of  differential  equations  make 
this  an  impractical  computational  method  for  solving  range- 
dependent  problems. 

In  an  attempt  to  reduce  the  complexity  involved  in  com¬ 
puting  the  coupling  matrices,  MeDonalds  used  the  original 
Picrcc-Mildcr  equations  to  argue  that  for  a  waveguide 
whose  horizontal  length  scales  arc  much  larger  than  the 


acoustic  wavelength  only  one  of  the  two  coupling  matrices 
has  significant  contribution.  By  neglecting  one  of  the  cou¬ 
pling  matrices  and  the  horizontal  derivative  of  the  density, 
McDonald  was  able  to  derive  an  expression  for  the  remain¬ 
ing  coupling  matrix  in  terms  of  local  modes  and  their  depth 
derivatives.  The  expression  for  the  coupling  matrix  derived 
by  McDonald  was  used  by  Abawi  et  al .6  who  derived  a  sys¬ 
tem  of  one-way  coupled  mode  equations  for  the  mode  am¬ 
plitudes  Although  this  method  is  a  practical  computational 
method  for  solving  range-dependent  problems,  it  still  suffers 
from  approximations  made  by  neglecting  one  of  the  two  cou¬ 
pling  matrices.  More  importantly,  since  this  method  is  based 
on  the  same  boundary  conditions  as  those  used  by  Pierce  and 
Milder,  the  energy  among  modes  is  not  conserved. 

The  coupled  mode  model  that  is  presented  in  this  paper 
is  derived,  not  from  the  wave  equation,  as  is  the  case  for  the 
Pierce  Milder  method,  but  from  the  equations  of  motion  for 
the  pressure  and  displacement  fields.  This  method,  which 
was  first  used  by  Shevchenko,  has  common  use  in  seismol¬ 
ogy  and  geophysics,  Odom,x  Maupin^  and  Tromp.10  The 
derivation  in  this  paper  follows  the  derivation  of  Tromp. 
While  the  model  derived  by  Tromp  is  for  the  general  clastic 
waveguide,  this  model  is  derived  for  a  waveguide  consisting 
of  fluid  layers.  Since  the  equations  of  motion  constitute  a 
system  of  two  first  order  coupled  differential  equations,  the 
coupled  mode  equations  resulting  from  them  only  contain  a 
single  coupling  matrix.  Furthermore,  this  method  provides  a 
natural  framework  for  applying  the  correct  boundary  and  in¬ 
terface  conditions  without  adding  any  more  complexity  to 
the  numerical  solution  of  the  equations.  In  fact,  it  is  shown  in 
this  paper  that  the  proper  application  of  the  boundary  and 
interface  conditions  not  only  simplifies  the  numerical  com¬ 
putation  of  the  coupling  matrix  by  allowing  it  to  be  ex¬ 
pressed  in  terms  the  local  modes  and  their  depth  derivatives, 
it  also  makes  it  possible  to  show  that  the  resulting  coupling 
matrix  is  anti-symmetric,  which  guarantees  the  conservation 
of  energy  among  modes.3,6,10 

The  method  presented  in  this  paper  and  the  references 
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cited  in  the  above  fall  in  the  category  of  continuous  eoupled 
mode  theory  where  the  solution  of  the  wave  equation  in  a 
range-dependent  waveguide  is  obtained  by  solving  of  a  set  of 
eoupled  differential  equations.  The  solution  for  a  range- 
dependent  waveguide  ean  also  be  obtained  by  the  discrete 
eoupled  mode  method.11  In  this  method  the  range-dependent 
waveguide  is  approximated  by  range-independent  stair  steps 
and  the  eoupled  mode  solution  is  obtained  by  matching  the 
solutions  of  the  wave  equation  for  neighboring  stair  steps  at 
their  common  boundary.  This  method  is  easy  to  implement 
numerically  and  has  wide  application  in  oeean  aeoustics. 
However,  the  method  that  is  presented  in  this  paper  has  two 
advantages  over  the  discrete  coupled  mode  method.  The  first 
advantage  is  it  expresses  the  coupling  matrix  in  terms  of 
physical  parameters  and  thus  provides  insight  into  the  pro¬ 
cess  of  mode  coupling  by  elearly  showing  what  is  respon¬ 
sible  for  it.  The  other,  more  important  advantage  is  that  this 
method  in  principle  ean  be  extended  to  handle  propagation  in 
three  dimensions,  where  the  discrete  eoupled  mode  method 
is  designed  for  propagation  in  two  dimensions  and  there  is 
no  obvious  way  to  modify  it  to  handle  propagation  in  three 
dimensions. 

This  paper  is  organized  in  the  following  way.  In  See.  II 
the  eoupled  mode  model  is  derived,  where  some  of  the  de¬ 
tails  of  the  derivation  are  given  in  the  appendices.  In  See.  Ill 
the  model  is  applied  to  two  range-dependent  propagation 
scenarios.  The  first  one  is  propagation  in  a  wedge,  where 
range-dependenee  is  entirely  due  to  variations  in  water 
depth.  The  parameters  used  in  this  example  are  sealed  to 
match  those  used  by  Coppens  and  Sanders1*  in  a  model  tank 
experiment.  The  seeond  example  is  propagation  through  in¬ 
ternal  waves,  where  range-dependenee  is  entirely  due  to  the 
variations  in  water  sound  speed.  The  parameters  used  in  this 
example  were  those  used  in  a  test  ease  in  the  Shallow  Water 
Aeoustics  Modeling  workshop.1  1  In  both  of  the  above  ex¬ 
amples  the  results  obtained  from  the  eoupled  mode  model 
are  compared  with  those  obtained  from  the  parabolic  equa¬ 
tion  (PE)  method.14 


II.  DERIVATION  OF  THE  COUPLED-MODE  EQUATIONS 

Consider  the  equations  of  motion  with  the  x-axis  in  the 
direction  of  propagation 

dxp  =  (uo2ux, 

dxltx— - 7  —  z' 

pc 

1 

z*u= - -  f)zp . 

pa) 

In  the  above  equations p  is  the  pressure  and  u  is  the  displace¬ 
ment  veetor.  The  pressure  and  the  normal  component  of  the 
displacement  are  continuous  aeross  any  interface.  This  may 
be  expressed  as 

[>]-=o,  [n-uj:=o. 


where  [£]!  =  £+  —  f  ,  the  +/—  indicate  the  value  of  the 
parameter  £  just  above/below  the  interface  and  n  is  the  unit 
veetor  normal  to  the  interface. 

The  field  quantities  ean  be  expressed  as  the  sum  of  nor¬ 
mal  modes 


=  S  P„{z)e'knX,  u(jr,e)=X  u„(x,z)eik"\ 

n  n 


where  pn  and  u„  denote  the  normal  modes.  Substituting  the 
above  expressions  into  the  equations  of  motion  results  in  the 
following  relationships  for  the  modes 


‘k„P„=po>2"„, 


Pn 


ikn»n= - - zd-_pn 


pc" 


pio 


1 

= - ~^zPn 


(i) 


bJ-=o,  [v«„]!=o. 

In  the  above  equations  un  denotes  the  x  component  of  the 
displacement. 

In  a  range-dependent  environment  the  pressure  and  the 
displacement  vector  may  be  expressed  as  a  sum  of  local  nor¬ 
mal  modes  with  range -dependent  coefficients,  cn(x) 


/>(*,z)  =  2  c„(x)pn(z-,x)ea"x, 

n 

(2) 

u(*>z)  =  2  c„{x)\in(:\x)e'k"'. 


In  this  notation  the  parametric  range-dependenee  of  the  loeal 
inodes  at  range  x  is  indicated  by  the  semieolon  separating  z 
and  x. 

Substitution  of  the  above  expansion  into  the  equations  of 
motion  gives 


cnPnelknX=^J  pC02CnU„elknX , 


Multiplying  the  first  equation  by  ume~lkn,x  and  the  seeond 
equation  by  pme~tknr\  adding  the  two  equations  and  inte¬ 
grating  along  the  depth  of  the  waveguide  gives 

2/  l  { i^nPnt  ~^~Pn^m^  ^x^n  ^  n(.P  m  ^x^n~^~^m  ^ xP  n ) 

n  JO 

+  icnkn{unpm+pnum)}e,{k”~k”')x  dz 
=  S  f - -p„Pm-Z-(d;U„)Pm 

n  JO  [  pc “ 

+  pa>2tinun^cne'(k”-k")*  dz.  (3) 

Sinee  according  to  Eq.  (I)  un= ik„p„/(p(o2),  we  find 
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2dx(cm)km+cm  <\km=  2  Amnc„. 


(4) 


The  details  of  the  above  derivation  is  given  in  Appendix  B. 
In  the  above  equation,  Amn  is  the  coupling  matrix  given  by 


ft*  1 

(km  +  k„)  ~Pm  <> xPn  dz 

Jo  P 

=  (k„-km)~]  Jo  |^|  dxi^  -jdj>„jp„ 


^  mn 


(kn  + 


CB\  CB 

^/w)  |  Pm  (.P kn  I  P nf 

Jo  P  Jo 


X<?jt|  -Jd~  +  l—P,tPmdxk 


,Hk„~km)x 


(5) 


The  above  equation  is  not  yet  in  the  desired  form,  as  it  con¬ 
tains  the  range  derivative  of  the  modes,  which  is  difficult  to 
compute  accurately.  In  the  remainder  of  this  section  we  will 
use  the  modal  equations  and  the  boundary  and  interface  con¬ 
ditions  to  convert  the  above  equation  into  one  which  only 
contains  the  local  modes  and  their  depth  derivatives. 

Consider  the  mode  equations  for  mode  n  and  mode  m 


+<?.  I j/jm+  -p„pmdx(k2) 


+  (*2-*'R[ -jP„Pm-^z\  -diPmjdtPnjdz-  (8) 

We  would  like  to  transfer  all  terms  involving  the  range  de¬ 
rivatives  of  the  modes  from  inside  the  integral  to  the  bound¬ 
ary  term  by  using  integration  by  parts.  The  first  term  can  be 
written  as 


JzPniPm 


<JS\  ~d:P»  I  4  (k2— k2)p„=0. 


+  -(k2-k2m)pm=0. 


(6) 

(7) 


The  modal  equation  is  obtained  by  substituting  un 
=  ik„pnJ(p<o2)  into  the  second  equation  in  Eq.  (1).  Next 
evaluate 


JM4) 

=  dx[-\(')zPn)Pm  ~f  dJ-}(dzPm)(d-P,Mz- 

L  \Pl  Jo  \Pl 

Similarly,  the  second  and  the  fifth  terms  can  be  written  as 
f”jz^z(^Pn)JP»,  dz 

(H  I 

'  ~9t  (p„)*:( 

Jo  P 


-di(9xp„)p„ 


zipJdz . 


and 


ru 

f*  /  1  \ 

{Pm  <?.»[Eq.  (6)]-[Eq.  (l)]PxPn}dz. 

Jo 

"  Jo  \  p**") 

This  gives 


~9x(Ph)(9zP* 


+ 


f»l 

~dz(Pn)d:(p„ 

Jo  P 


)dz. 


dzPn  )Pm  +  f>z\  ~dz(dxpa)  \pn 


+  -PnPm  dx(k2)  +  {k2~  k2„)r)x\  ~)PnPn 


1 


Substituting  these  into  b'q.  (8)  gives 

f»  I 

(km+kn)  -p„,  dxp„dz 
Jo  P 


+  l k;„-kn)-p„,  r')xPn-  -  d-_pm  )  dxp„ 


dz  —  0. 

The  fifth  term  in  the  above  equation  can  be  written  as 


+  ~PnPm  dx(k2)  +  (A-2  - -jp„pm  )  dz 


[» 1 

^«)( kn)  I  Pm  ^ xP n  *^Z 
Jo  P 

=  “  Jn  |  dzPn  )p,„+  “  <?i( djrPn )  jp» 


-J(d.p„)p„+  -dx{d-p„)pn 


■dx(pn)(d.pm) 


(9) 


+  -PnPm  d.x (k2)  +  ( k 2  -  -JPnPn 


dz. 


The  boundary  term  in  the  above  equation  may  be  written  as 


-Oz\-PzP,nj<hp„ 


<M  -J(JzPn)PmJ-  -dx(dspn)pm-  -dx(pn){azpm) 


d.x\  ~°-Pn  \Pm-dx(p„)-((hpJ 
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l(k2-k„km)p„pm 


Sinee  the  derivative  along  the  interface  of  a  continuous  func¬ 
tion /is  continuous,  we  have 


(kni  ktl)Amn 


*  *  dh„ 

[T  V/]  =0,  where  T—  \+  — z. 


^zPn  dzPm)dX\  ~)  +  x(k2) 

Pi  P 


dz 


This  gives 

[**/!-  = -[**(*)**(/)]-■ 

Since  both  pn  and  dzp„/p  are  continuous  across  the  interface 
z  =  h(x),  the  above  boundary  term  can  be  written  as 

1  ll 

dxWd.(pn)~(dj>m)  ~  dxh  d.y-dtpn 

With  the  help  of  the  wave  equation,  Eq  (6),  this  may  be 
written  as 


dx(h)a2(p„)-(d:Pm)  +  dxh-(k2-kl)p„p„ 


Substituting  this  into  Eq.  (9)  and  the  result  into  Eq.  (5)  yields 


4-  drh 


-dzPn  dzpm+~{k2 


~k„km)pnp„ 


i  (10) 


The  expression  for  the  coupling  matrix  given  by  the  above 
equation  is  the  main  result  of  this  paper.  It  shows  the  effeet 
of  mode  coupling  due  to  contribution  from  volumetrie  and 
bathymetrie  variations  in  the  waveguide  separately.  The  first 
part  of  the  coupling  matrix  containing  the  integral  is  due  to 
contribution  from  volumetrie  variations  in  the  waveguide 
such  as  variations  in  sound  speed  and  density  as  a  function  of 
range.  The  seeond  part  is  due  to  contribution  from  bathymet¬ 
rie  variations  in  range,  i.e  ,  variations  in  water  depth,  as  is 
evident  from  the  presence  of  dxh. 

The  coupling  matrix  has  two  important  properties.  First, 
it  is  anti-symmetrie,  i.e.,  Am„=  —  A  .  This  implies  that  en- 


Coupled  Mode  Solution 


Range  (km) 


Adiabatic  Mode  Solution 


Depth  Slice  at  30  m 


Range  (km) 


FIG.  I.  Propagation  in  an  oceanic  wedge:  The  water  depth  is  constant  at  200  in  for  the  first  5  km  and  it  slowly  decreases  tu  zero  in  the  next  7.5  km  resulting 
in  a  wedge  angle  of  1.55  deg.  The  top  left  panel  shows  the  acoustic  field  due  a  25  Hz  source  placed  at  180  m  as  a  function  range  and  depth  computed  using 
the  coupled  mode  model.  The  bottom  left  panel  shows  the  acoustic  field  computed  using  the  adiabatic  mode  model.  The  top  and  bottom  panels  on  the  right 
show  a  comparison  of  the  transmission  loss  as  a  function  of  range  computed  using  the  PF  model  (solid)  the  coupled  mode  model  (dotted)  for  two  receiver 
depths. 
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FIG.  2.  Propagation  through  internal  waves:  The  ocean  environment  consists  of  a  200  m  layer  over  a  bottom  half  space.  The  sound  speed  in  the  water  is 
modeled  to  simulate  fluctuations  due  to  internal  waves.  The  top  left  panel  shows  the  sound  speed  fluctuations  in  the  water.  The  top  right  panel  shows  the 
acoustic  field  in  this  environment  due  to  a  100  Hz  source  placed  at  100  ni  computed  using  the  coupled  mode  model.  The  bottom  two  panels  show  a  comparison 
of  the  transmission  between  the  PE  model  (solid),  the  coupled  mode  model  (dotted),  and  the  adiabatic  coupled  mode  model  (dashed).  The  transmission  loss 
is  computed  at  a  receiver  depth  of  100  m  and  for  source  frequencies  of  25  and  100  Hz. 


ergy  is  conserved  among  modes.  Second,  it  only  contains  the 
modes  and  their  depth  derivatives.  This  means  that  the  cou¬ 
pling  matrix  can  easily  be  computed  using  the  local  modes 
and  their  depth  derivatives,  which  can  be  obtained  from  any 
normal  mode  code  such  as  KRAKEN. ls 

111.  EXAMPLES 

In  this  section  the  one-way  coupled  mode  model  devel¬ 
oped  in  the  previous  section  is  applied  to  two  range- 
dependent  cases.  In  the  first  example  we  use  the  one-way 
coupled  model  to  compute  acoustic  propagation  an  oceanic 
wedge  wherein  range-dependence  is  due  to  variations  in  the 
water  depth.  In  the  second  example  propagation  through  an 
ocean  with  internal  waves  is  computed  where  the  ocean  en¬ 
vironment  is  chosen  such  that  rangc-dcpcndcncc  is  entirely 
due  to  variations  in  sound  speed.  The  results  arc  compared 
with  those  obtained  using  the  parabolic  equation  PE14 
method. 

A.  Propagation  in  a  wedge 

The  ocean  environment  in  this  example  is  scaled  to  cor¬ 
respond  to  the  model  tank  experiment  reported  by  Coppcns 
and  Sanders.12  The  water  depth  is  initially  200  m  for  the  first 
5  km  and  then  it  slowly  decreases  to  zero  in  the  next  7.5  km 
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resulting  in  a  wedge  angle  of  approximately  1.55  deg.  To 
approximate  the  branch  cut  integral  in  the  modal  representa¬ 
tion  of  the  field,  a  1000  m  deep  false  bottom  is  used.  The 
waveguide  consists  of  two  isovelocity  layers:  a  water  layer 
over  a  bottom  layer.  The  water  sound  speed  is  1500  m/s  and 
its  density  is  1  0  g/cm"  The  bottom  sound  speed  is  1700  m/s 
with  a  density  of  1.15  g/cm3.  The  attenuation  in  the  bottom 
is  0.5  dB/\.  A  25  Hz  source  is  placed  at  180  m.  These 
environmental  parameters  arc  chosen  to  correspond  to  those 
used  by  Coppcns  and  Sanders. 

The  acoustic  field  in  the  waveguide  is  computed  using 
Eq.  (2)  with  the  modal  coefficients,  cm ,  obtained  from  the 
solution  of  Eq.  (4).  The  first  order  differential  equation  for 
the  modal  coefficients,  Eq.  (4),  is  solved  by  using  fourth- 
order  Rungc  Kutta  integration  To  obtain  the  modes  and  the 
coupling  matrix  as  a  function  of  range,  the  wedge  is  divided 
into  range-independent  stair  steps.  The  local  modes  and  the 
local  coupling  matrix  using  Eq.  (10)  arc  computed  in  each 
stair  step  and  updated  in  the  differential  equation  The  step 
size  in  the  two  examples  that  arc  presented  in  this  paper  is  10 
m  However,  a  step  size  of  50  m  gives  identical  results. 

The  results  of  the  above  computation  are  shown  in  Fig. 
1.  The  top  left  panel  in  Fig.  1  shows  the  acoustic  field  com¬ 
puted  using  the  one-way  coupled  mode  model  described  in 
this  paper.  It  can  be  seen  that  as  the  water  depth  decreases, 
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the  water  modes  (there  are  three  water  modes  in  this  ex¬ 
ample)  cutoff  in  the  form  of  discrete  beams  radiating  into  the 
bottom.  The  experimental  data  obtained  by  Coppens  and 
Sanders  show  the  exact  same  behavior.  Jensen  and 
Kuperman1*  used  the  parabolic  equation  method  to  model 
the  acoustic  propagation  in  this  example  and  found  results 
identical  to  those  shown  in  the  top  left  panel  of  Fig.  1.  They 
interpreted  the  slow  disappearance  of  the  discrete  water 
modes  into  the  bottom  as  an  indication  that  energy  contained 
in  a  given  mode  docs  not  couple  into  the  next  lower  mode 
but  couples  almost  entirely  into  the  continuous  mode  spec¬ 
trum.  While  this  effect,  which  is  a  consequence  of  mode 
coupling,  is  implicitly  accounted  for  in  the  parabolic  equa¬ 
tion  (PE)  formulation,  the  coupled  mode  model  explicitly 
accounts  for  it  through  the  coupling  matrix,  The  two  panels 
on  the  right  in  Fig.  1  show  a  comparison  of  the  transmission 
loss  computed  using  the  one-way  coupled  model  and  the  PE 
model  for  two  receiver  depths.  The  close  agreement  between 
the  two  models  clearly  demonstrates  that  the  one-way 
coupled  mode  model  correctly  accounts  for  mode  coupling. 
If  the  coupling  matrix  in  the  one-way  coupled  mode  model 
are  set  equal  to  zero,  the  one-way  coupled  mode  model  re¬ 
duces  to  the  adiabatic  mode  model.  The  bottom  left  panel  in 
Fig.  1  shows  the  acoustic  field  computed  using  the  adiabatic 
mode  solution.  Observe  that  the  adiabatic  mode  solution 
docs  not  have  the  correct  field  behavior  near  cutoff.  While  in 
the  coupled  mode  solution  modes  gradually  radiate  their  en¬ 
ergy  to  the  bottom  near  cutoff,  in  the  adiabatic  mode  solution 
this  process  occurs  abruptly  because  there  is  no  mechanism 
for  the  transfer  of  energy  between  modes. 


B.  Propagation  in  internal  waves 

The  ocean  environment  in  this  example  is  one  of  the  test 
eases  used  at  the  Shallow  Water  Acoustic  workshop.1 5  It  con¬ 
sists  of  200  m  of  water  over  a  400  m,  isovcloeity  bottom. 
The  bottom  density  1.5  g/cm3  and  its  sound  speed  was  1700 
m/s.  The  sound  speed  profile  and  the  velocity  fluctuations 
due  to  internal  waves  in  the  water  column  arc  modeled  ac¬ 
cording  ton 

c(:/)  =  c(z)  +  Ac(z,r), 
where 


c(z)  = 


151 5.0+0.01 6z  z<26 
1490(  1 .0+  0.25(e  h  +  b  —  1 .0)) 


z>26 


and 


A c:(z,r)  =  —e  z/25cos  27rr. 

In  the  above  b  =  (z  —  200)/500  and  r  is  measured  in  km.  The 
top  left  panel  in  Fig.  2  shows  the  sound  speed  profile  for  this 
example.  The  top  right  panel  shows  the  acoustic  field  in  the 
waveguide  for  a  100  FIz  source  place  at  30  m  The  bottom 
two  panels  in  Fig.  2  show  a  comparison  of  the  adiabatic 
normal  mode  and  the  coupled  mode  solutions  with  the  PE 
solution  for  two  source  frequencies.  In  both  eases  the  re¬ 
ceiver  depth  is  at  70  m.  The  coupled  mode  solution  agrees 
well  with  PE  solution  at  both  frequencies  while  the  adiabatic 


normal  mode  solution  docs  not  agree  with  the  parabolic 
equation  solution  at  all.  This  is  more  evident  at  the  higher 
frequency  where  the  effects  of  the  internal  waves,  and  thus 
the  mode  coupling  due  to  them,  is  stronger. 
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APPENDIX  A:  THE  MODE  ORTHOGONALITY 
RELATIONSHIP 


The  mode  orthogonality  equation  is  obtained  by  multi¬ 
plying  the  mode  equation  (6)  by  pm  and  Eq.  (7)  by  pn  sub¬ 
tracting  the  resulting  equations  and  integrating  to  give 


dz 


dz  —  0. 


Integrating  the  first  integral  by  parts  results  in  boundary 
terms 


B 


Pm-dzPn 


-Pn-^zPn, 
0  * 


0 


Since  either  the  mode  or  its  derivative  is  zero  at  the  bound¬ 
aries,  there  is  no  contribution  from  the  above  interface  terms. 
What  remains  is 


f*  1 

(km-kn){km  +  kn)  -pmp„d:  =  Q. 

Jo  (> 

For  m  =£  /?, 

t*  1 

(km+k„)\  -pmpndz  =  0, 

Jo  P 

and  for  m  —  n  we  choose  to  normalize  the  inodes  such  that 

[B  1 

(*« + *„)  ~pmp „dz  =  2k„ S„m . 

Jo  P 


APPENDIX  B:  DETAILS  LEADING  TO  EQ.  (4) 

We  start  by  multiplying  the  first  equation  in  Eq.  (1)  by 
um  and  the  second  equation  in  Eq.  (1)  by  pm  ,  adding  the  two 
equations  and  integrating  to  get 

‘k„  J  (pnum  +  unpm)dz 
Jo 


PnPm 
pc 2 


1 

- ~,flzPn 

f)CO~ 


Pm+po>-unu„,\d:. 
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Integrating  the  middle  term  on  the  right-hand  side  by  parts 
yields 

\ 


Jo  \pc 2  / 


Pm  d-=\  I  — 

poj 


-n-w 

Jo  \  p(0* 


d.pm  dz. 


Since  the  quantities  inside  the  square  brackets  are  continuous 
across  the  interface,  the  boundary  term  is  zero.  This  results  in 


2  +  - W 


,  (B,  PnPm 

‘K„  (Pnum  +  U„Pm)dz= - - 

Jo  Jo  \  pc 


pc  phi 


X(V>m)  +  PW2U„U„ 


dz.  (Bl) 


Next  consider  the  following  term  in  Eq.  (I),  which  can  be 
integrated  to  give. 


f  Z'(V„ 

Jo 


)Pm  clz  =  [  ™nPmY- 


I  (z-u«  da 

Jo 


Pm)dZ- 


(B2) 

The  continuity  condition  for  the  normal  component  of  the 
displacement  can  be  written  as 

[n-u„]^  =  [(  — z+  djix)  •  u„]*  =  0. 

This  gives, 

[z-u 

Since  pm  is  continuous  across  the  interface,  Eq.  (B2)  reduces 
to 

f  i-(dzun)pm  dz=[(dxh)unpm]t~  f  (z-u nf)zPm)dz. 

Jo  Jo 

Substituting  this  into  Eq.  (3)  and  using  Eq.  (BI)  gives 

I  ^XC n  f  ( u np p nit m)dz  -f  c n  I  (/? m  0xun 
n  [  JO  Jo 

Next  substituting  for  un=  iknpnl  urp  gives 

2  (dxc„)(k„  +  km)  f  -p„pm  dz+c„(k„+km) 

n  JO  P 

Pm  fB  1 

^  I  (  ^xP n)  dz  ^~Cn{  dxkn)  I  P nP ru  dz 

Jo  P  Jo  p 

Jnk"j0  ~|  PnPmdz 


+  cn 


-c„\  (dxh) 


knPnPn 


ei(kn-km).x\=0 


(B3) 


Using  the  orthogonality  of  the  modes  we  find, 
2  c m  ^x^m~  S  ^mncni 


where  the  coupling  matrix,  /lww  is  defined  by  Eq.  (10).  For 
m  —  n  Eq.  (B3)  becomes. 


2  dxcnkn+cn  dxk 


.  +  cjt.j0'(  2  (rJxPn)j+p2„  dx\  ~l  I  dz 


(VO 


knP'„ 


=  0. 


This  can  be  written  as 


f  ^  i  P” 

^  n  ^xkn~\~  Cnkn  I  dx  I 

Jo  \  P 


+  C. 


{dxh) 


KpI 


■  0. 


It  is  shown  in  Appendix  C  that 

•  fiYh 


which  gives 

2  dxc„k„+c„  dxk„  =  0. 


—pi 

P,n 


(B4) 


APPENDIX  C:  DERIVATION  OF  EQ.  (B4) 

The  integral  across  the  waveguide  can  be  written  as, 

/  »2\ 


IMvh-M  jH'Ajb 


Each  one  of  the  above  integrals  can  be  written  as 

,,2\  r.,l„  2^ 

V) 


and 


b  j  p2 

cli,n 


Substituting  for  these  quantities  we  find 


The  integral  on  the  right-hand  side  is  a  constant  which  gives 
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